Quasiparticle transport in the vortex state of YBa 2 Cu 3 06.9 
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The effect of vortices on quasiparticle transport in cuprate superconductors was investigated by 
measuring the low temperature thermal conductivity of YBa2Cu30e.9 in magnetic fields up to 8 T. 
The residual linear term (as T — > 0) is found to increase with field, directly reflecting the occupation 
of extended quasiparticle states. A study for different Zn impurity concentrations reveals a good 
agreement with recent calculations for a d-wave superconductor. It also provides a quantitative 
measure of the gap near the nodes. 

PACS numbers: 74.70.Tx, 74.25.Fy 

In unconventional superconductors, where the gap 
structure has nodes along certain directions, it was ar- 
gued by Volovik H that the dominant effect of a magnetic 
field should be the Doppler shift of extended quasiparti- 
cle states due to the presence of a superfiuid flow around 
each vortex. Reports of a VH field-dependence in the 
heat capacity of YBa 2 Cu306.9 (YBCO) have generally 
been accepted as strong evidence for this effect [|| . How- 
ever, a similar field dependence is also observed in an 
s-wave superconductor such as NbSe2 ||, where it can 
only arise from localized states bound to the vortex core. 

In this respect, it is interesting to look at heat con- 
duction, to which only delocalized states contribute. So 
far, in all measurements performed above 5 K or so, the 
thermal conductivity k is found to drop with field [^-^), 
much as it does in a clean type-II superconductor like 
Nb 0, and eventually levels off to a roughly constant 
plateau In Bi 2 Sr 2 CaCu20 8 (BSCCO), the crossover 
from drop to plateau occurs abruptly, at a field whose 
magnitude increases with temperature Q . Franz has re- 
cently proposed that the plateau is due to a compensation 
between the increasing occupation of extended quasipar- 
ticle states a la Volovik and a parallel increase in the 
scattering of quasiparticles by vortices 0. While this 
may well be part of the explanation, it can hardly serve 
as a direct confirmation of the "Volovik effect" . Neither 
can it explain the abruptness of the onset of the plateau. 
The situation remains quite puzzling, in part because of 
the large phonon background at these temperatures. 

In this Letter, we present a study of heat trans- 
port which provides a solid experimental basis for a 
description of quasiparticle properties in terms of the 
field-induced Doppler shift of a d-wave spectrum due to 
the superflow around vortices. Furthermore, the good 
agreement found with recent calculations by Kiibert and 
Hirschfeld |[o| allows us to conclude that in YBCO at 
low temperature, impurity scattering is close to the uni- 
tarity limit and vortex scattering is weak. Moreover, it 
provides a measure of the two parameters that govern 



the Dirac-like spectrum near the nodes, responsible for 
all low-energy properties. 

We have measured the thermal conductivity of 
three optimally-doped untwinned single crystals of 
YBa2(Cui_ a; Zn a ;)306.9 down to 0.07 K, with a current 
along the a-axis of the basal plane and a magnetic field 
along the c-axis. The nominal Zn concentrations are 
x = (pure), 0.006 and 0.03, corresponding to T c = 93.6, 
89.2, and 74.6 K, respectively. The data were taken by 
sweeping the temperature at fixed field, and the field was 
changed at low temperature (below 1 K) . Data taken by 
field-cooling the samples gave the same results. Given 
that the field of full penetration is estimated to be about 
1 T the lowest field point was set at 2 T. The largest 
source of absolute error is the uncertainty in the geomet- 
ric factor (approximately ± 10%); however, the relative 
error between different field runs on the same crystal is 
only due to the fitting procedure used to extract the elec- 
tronic contribution (see below), which is less than 10%. 

Fig. 1 shows the measured thermal conductivity di- 
vided by temperature, k/T, as a function of T 2 in fields 
up to 8 T. As discussed in Ref. (l2), the only reliable way 
to separate the electronic contribution from the phonon 
background is to reach the T 3 regime for the phonon 
conductivity (below about 130 mK), from which a linear 
extrapolation to zero yields the purely electronic term at 
T=0, labelled k„/T. 

Zero field. — Let us first focus on the results at H=0. A 
linear term is observed at T=0, of comparable magnitude 
in all three crystals. From measurements on a number of 
pure (optimally-doped) crystals of YBCO, we obtain an 
average value for J 1 1 a Il3) : 



T 



0.14 ±0.03 mW K" 



(1) 



This residual conduction is associated with an impurity 
band whose width 7 is a new energy scale relevant to all 
low-energy properties in superconductors with gap zeros 
p3 |. 7 grows with the impurity scattering rate T, in a 
way which depends strongly on whether impurities act 
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FIG. 1. Thermal conductivity divided by temperature vs 
T 2 of YBa2(Cui_ a; Zn a ;)306.9 in a magnetic field, for H = 
(triangles), 2 (diamonds), 4 (circles), and 8 T (squares). 
The bottom panel corresponds to x~0; middle, £=0.006; top, 
£=0.03. The lines are linear fits to the data below 130 mK. 

as Born or resonant scatterers. In the unitarity limit, 
7 is maximum and given by 7 ~ 0.61\/TlTAq, where 
Ao is the gap maximum p9[ . This novel fluid, deep in 
the superconducting state, should carry heat and charge. 
Lee was the first to point out the universal character of 
that conduction as T — > (and frequency u> — > 0), 
for certain gap topologies Pq |. For the d x 2_ y 2 gap, the 
conductance of a single CuC>2 plane is predicted to be 
coo = ^pf- \ where vp and v 2 are the two pa- 
rameters governing the Dirac-like spectrum near a node: 
E(k) = y/e{ + A 2 k = hy/v'jrki + v'm. Here (fci.fa) de- 
fines a coordinate system whose origin is at the node, 
with k\ normal to the Fermi surface and k 2 tangen- 
tial; thus vp is the Fermi velocity and v 2 is the slope 
of the gap, at the node. The conductivity is independent 
of impurity concentration as a consequence of the ex- 



act compensation between the decreasing scattering time 
t oc I/7 and the growing residual quasiparticle density 
of states, N(0) oc j/(vpv 2 ), when impurities are added, 
i.e. coo « N(0)v f t cx vp/v 2 . 

Graf and co-workers showed that the residual normal 
fluid should obey the Wiedemann- Franz law at T=0 Q , 
so that the thermal conductivity as T — ► in a rf-wave 
superconductor should also be universal, and given by: 



k'i 



Vp 



= L (Too = tt| ~ n ' 
1 6n v 2 



(2) 



where n is the number of CuC>2 planes per meter stacked 
along the c-axis and Lq = \{^f-) 2 - 

In recent measurements of n(T) in YBCO, the univer- 
sal limit was observed fl2]| . It is also evident from the 
H=Q curves in Fig. 1, noting that the scattering rate T 
in the 3% Zn crystal is some 20 to 40 times larger than it 
is in the pure crystal. Indeed, estimates from a combina- 
tion of resistivity, microwave and infrared measurements 
give hT/k B T c0 = 0.014, 0.13, and 0.54 for x = 0, 0.006 
and 0.03, respectively Note also that the impurity 
bandwidth in the 3% Zn sample is a sizable fraction of 
the gap maximum, so that corrections to the universal 
limit are expected. Calculations by Sun and Maki give a 
30% increase in kq/T for 20% T c suppression (ll| , in good 
agreement with the observed slight increase (see Fig. 2). 

For the bi-layer structure of YBCO, n — 2/c where 
c = 11.7 A is the c-axis lattice constant. Combining this 
with Eqs. (1) and (2) yields 

— = 14±3. (3) 

Taking Vp = 1.2xl0 7 cm/s (see Ref. 0), one gets 
v 2 — 1 x 10 6 cm/s, which corresponds to a slope of the gap 
at the node dA((j>) /d<j> = hkpv 2 — 380 K, very close to 
that estimated using the simplest d-wave gap, Aocos2</>, 
in the weak-coupling limit: 2A = 2 x 2.lAk B T c = 400 K. 
The same ratio vp/v 2 governs the linear drop in super- 
fluid density with temperature when 7 <C k B T <C Ao. In 
terms of penetration depth data, where A(T) = A(0) + 
5X(T), one has SX(T) = (^)2 In 2(^) @, where 



2ir 2 kpv 2 <T 0( ), so that 8X/T - 



4.8 A/K for the above 



values of vf and v 2 . This is in excellent agreement with 
the microwave measurements of Zhang and co-workers, 
who obtain SX a /T = 4.7 A/K 

Vortex state. — As seen in Fig. 1, the effect of a mag- 
netic field applied normal to the Cu02 planes is to in- 
crease K a in YBCO. To ensure that only the electronic 
contribution is considered, the T — ► limit of k/T is 
plotted as a function of field in Fig. 2, for each of the 
three crystals. The enhancement observed as T — > is 
direct evidence for itinerant quasiparticle excitations in- 
duced by the field. This is in sharp contrast with the 
behaviour found both in conventional superconductors 
M and in YBCO at higher temperatures where 
low fields (compared to H C 2) suppress the conductivity. 



2 




H(T) 

FIG. 2. Residual linear term Ko(H)/T normalized by 
Koo/T as a function of applied field for pure (squares), 0.6% 
Zn (circles) and 3% Zn (triangles) samples. Fits to Eq.(4) for 
each crystal yield the values of p shown (in units of T~ 1,/2 ). 

Volovik has shown how even at T=0 a small mag- 
netic field can produce quasiparticles near the nodes 
of a cZ-wave gap In the presence of a superfluid 

flow, the quasiparticle spectrum is Doppler shifted to 
E(k, A) = E(k) — |vk ■ A, where Vk is the normal state 
velocity and A is the vector potential jl7) . In the case of 
vortices, the magnitude of the Doppler shift may be char- 
acterized by its average, Eh, obtained by integrating over 
a vortex-lattice unit cell (of radius R) : Eh = f (vk • A) ~ 
■jp J dr r vp-A ~ ^$>q>vf ~ VH, using 2R ~ intervortex 
spacing ~ 1/ yH and the fact that <I>o = £ A • d£, with $o 
the flux quantum. Volovik first derived the y/~H depen- 
dence and obtained the following field-induced density 
of states (DOS) for a <i-wave gap in the clean limit [Q: 
SN(0;H)/N ~ y/8/TT E H /A , where N is the DOS at 
the Fermi level, Eh — ah 



vp y/H/^o and a is a 
vortex-lattice parameter of order unity jl9| . The effect is 
sizable even at low fields; in YBCO at 8 T, one expects 
E H ^ 40 K and SN(0; H)/N„ ~ 30% (for o=l). 

Comparison with theory. — Kubert and Hirschfeld have 
calculated the electronic thermal conductivity of a d-wave 
superconductor in a field, neglecting vortex scattering 
pi. At T=0, they obtain for H 1 1 c and J _L c 



k(0;H) 
f 



T /VI 



sinh 



(4) 



where p is essentially the ratio of the two relevant en- 
ergy scales, 7 and Eh\ in the dirty limit where Eh < 7, 
p = \[§pK 7 1 Eh oc 1/^/H. A one-parameter fit of the 
data to Eq.(4) for each crystal is shown in Fig. 2. The 
first point to note is that the sub-linear dependence on 
field is well reproduced. Perhaps more important is the 
fact that the magnitude of the response in all three cases 
is very much as expected. Indeed, the fits yield the fol- 
lowing values for p, evaluated at 8 T: 0.92, 1.63, and 
4.32 for x = 0, 0.006, and 0.03, respectively. These cor- 
respond to a ratio j/Eh = 0.67, 1.18 and 3.12, which 
shows that none of the crystals is in the clean limit over 
the field range investigated. Treating impurity scatter- 
ing in the unitarity limit, the scattering rate becomes: 
hT/k B T c0 = 0.02, 0.07, and 0.5, respectively (taking 
A = 2.Uk B T c0 ), given that E H ~ 20 K at 8 T (as- 
suming vp = lxlO 7 cm/s and a=l/2). These values are 
in remarkable agreement with those given above as inde- 
pendent estimates based on the residual resistivity. 

The same theoretical treatment has been applied to 
the specific heat. Kubert and Hirschfeld have obtained 
A (()://;- A i, ~ 



^7*/7n = y / ^/ 7ra V H / H c2, where 5-f* is 
the increase in specific heat due to an applied field and 
7 n is the normal state specific heat at the Fermi level 
p9| . We can re- formulate this expression in terms of a 

and v 2 such that Sj* = pVH where f3 = ^f--^^. 
So now we have three equations for the three parame- 
ters a, vp and v 2 - Using (3 = 0.9 mJ K~ 2 T -1 / 2 mol -1 
§j, we find v 2 = 2.2axl0 6 cm/s. From k o (0;H)/T, we 
have avp = 5xl0 6 cm/s. These two expressions, to- 
gether with Eq.(3), yield a = 0.40. With this value of a, 
v 2 = 0.9xl0 6 cm/s and vp = 1.2xl0 7 cm/s. All three 
values are very reasonable, so that the overall picture is 
quantitatively convincing, lending strong support to the 
basic mechanism of a Doppler shift by superflow around 
vortices. More generally, it validates the theory of trans- 
port in unconventional superconductors; in particular, 
the assumption that in these correlated electron systems 
scattering by impurities must be treated in the unitarity 
limit of strong, resonant scattering appears to be correct. 

Other measurements. — Let us now compare our re- 
sults with other k measurements. Note first that previ- 
ous low temperature measurements in YBCO were done 
with the field in the basal plane where the Volovik effect 
is expected to be much smaller. The two studies were 
in disagreement, with a 50% increase in 8 T in one case 
Ipof , but no change detected in 6 T in the other ^3] . At 
higher temperatures, the fact that k decreases with field 
could come from a number of effects. Perhaps the most 
natural is vortex scattering, as invoked in the case of Nb. 
While this can apply to both quasiparticles and phonons, 
the fact that the largely electronic peak below T c is al- 
most completely suppressed in 10 T [Q suggests that the 
former suffer most of the impact. Franz has shown that 
a disordered vortex lattice in a d-w&ve superconductor 
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can result in quasiparticles scattering off the superflow 
P|. At low temperature, however, we find no indication 
of significant vortex scattering, in the sense that there is 
good agreement between calculations neglecting vortex 
scattering and our data on crystals for which the relative 
strength of impurity and vortex scattering must differ 
markedly from one crystal to the next. Further work is 
needed to arrive at a coherent description of quasiparticle 
transport in the vortex state which covers both regimes 
of behaviour. 

It is interesting to compare our results on YBCO with 
the corresponding results on BSCCO, obtained recently 
by Aubin, Behnia and co-workers 0]. At tempera- 
tures below 0.7 K, an increase in k/T with field is also 
found, again with a sublinear (roughly \f~H) dependence. 
What is striking is the magnitude of the response. The 
application of only 2 T increases k/T by about 0.20 



mW K 2 cm 



Now from our own measurements on 



pure, optimally-doped single crystals of BSCCO, the 
(presumably universal) residual linear term is kq/T = 
0.18±0.03 mW K" 2 cm" 1 This means that in 

BSCCO a field of 2 T causes the quasiparticle conduction 
to double, whereas it produces only a 35% increase in our 
pure YBCO crystal. In reality, the field dependence in 
BSCCO is much more than 3 times stronger, since the 
impurity scattering rate in the crystal used by Aubin and 
co-workers could be as much as 100 times larger. Indeed, 
its residual resistivity is 130 fj.fl cm p3], compared with 
approximately 1 /ifi cm in our pure YBCO crystals Jl2| . 
These considerations lead us to conclude that the nature 
of defect scattering in these two (otherwise quite simi- 
lar) materials is significantly different. Either the kind 
of defect found in nominally pure crystals is different, or 
the impact that a given defect (e.g. impurity) has on the 
surrounding electron fluid is different. It will be inter- 
esting to see whether the current theoretical framework 
can account for the BSCCO data simply by adjusting the 
impurity phase shift, moving it away from the unitarity 
limit towards the Born limit. 

In summary, we have measured the low temperature 
heat conduction by quasiparticles in the Cu02 planes of 
YBCO as a function of magnetic field, for different im- 
purity levels. In all cases, the residual linear term kq/T 
is found to increase with field strength, directly reflect- 
ing the additional population of extended quasiparticle 
states. The good agreement with calculations by Kiibert 
and Hirschfeld for a d-wave superconductor allows us to 
draw the following conclusions: 1) the "Volovik effect" is 
fully verified, and it is the dominant mechanism behind 
the field dependence of transport in YBCO as T — > 0, 
for H || c; 2) vortex scattering, invoked to explain the 
behaviour at intermediate temperatures, is seemingly ab- 
sent at low temperature; 3) the widespread assumption 
that impurities (or defects) can be treated as unitary 
scatterers in correlated electron systems appears to be 
verified in YBCO; 4) the nature of impurity scattering 



differs between BSCCO and YBCO. 

We thank K. Behnia, P. Hirschfeld and P. A. Lee for in- 
sightful discussions, and B. Lussier for his help during the 
early stages. This work was supported by the Canadian 
Institute for Advanced Research, and funded by NSERC 
of Canada and FCAR of Quebec. MC acknowledges sup- 
port from the Carl Reinhardt Foundation. 

t Current address: Department of Physics, University of 
Toronto, Toronto, Canada M5S 1A7. 



[3 

[4 
[5. 

[6] 

[r] 

[8. 
[9 
[10 

[11 
[12 
[13 



[14 
[15 
[16 
[17 



[18 
[19 

[20 
[21 
[22 
[23 
[24 



G. E. Volovik, J. E. T. P. Lett. 58, 469 (1993). 
Kathryn A. Moler et al, Phys. Rev. Lett. 73, 2744 

(1994) ; ibid, Phys. Rev. B 55, 3954 (1997); D. A. Wright 
et al., Phys. Rev. Lett. 82, 1550 (1999). 

D. Sanchez et al Physica B 204, 16 7 (1995); J. E. Sonier 
et al., preprint ( x>nd-mat/981142C ). 
T. T. M. Palstra et al., Phys. Rev. Lett. 64, 3090 (1990). 
S. D. Peacor, J. L. Cohn and C. Uher, Phys. Rev. B 43, 
8721 (1991). 

K. Krishana et al., Science 277, 83 (1997). 

J. Lowell and J. B. Sousa, J. Low Temp. Phys. 3, 65 

(1970). 

K. Krishana et al., preprint (1998). 

M. Franz, Phys. Rev. Lett 82, 1760 (1999). 

C. Kiibert and P. J. Hirschfeld, Phys. Rev. Lett. 80, 4963 
(1998). 

Ernst Helmut Brandt, Physica C 282-287, 343 (1997). 
Louis Taillefer et al., Phys. Rev. Lett. 79, 483 (1997). 
This average of three measurements, taken with a greater 
point density at low temperature, is more accurate than 
the previously reported 0.19L0.03 mWK" 2 cm _1 [12]. 
(Note that the two values do agree within error.) 
See, e.g., M. J. Graf et al., Phys. Rev. B 53, 15147 (1996). 
P. A. Lee, Phys. Rev. Lett. 71, 1887 (1993). 
Y. Sun and K. Maki, Europhys. Lett. 32, 355 (1995). 

D. Xu, S. K. Yip and J. A. Sauls, Phys. Rev. B 51, 16233 

(1995) ; P. A. Lee and X.-G. Wen, Phys. Rev. Lett. 78, 
4111 (1997). 

, Phys. Rev. Lett. 73, 2484 (1994). 
J. Hirschfeld, Sol. St. Commun. 105, 



Kuan Zhang et al. 
C. Kiibert and P. 
459 (1998). 
CD. Bredl et al., 
B. Wand et al., J. 
Herve Aubin et al 



Z. Phys. B 86, 187 (1992). 
Low Temp. Phys. 105, 993 (1996). 
, Phys. Rev. Lett. 82, 624 (1999). 
Patrik Lambert, M.Sc. Thesis (McGill University, 1998). 
K. Behnia, private communication. 



4 



